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ABSTRACT 

In this letter, we investigate the nature of the dramatic 'super-bubble' emission struc- 
tures associated with the z=2.44 radio galaxy MRC 0406-244, using rest-frame optical 
spectroscopy and an archival HST NICMOS image. Based on the optical line ratios 
and the HST morphology, we conclude that the gas in the superbubbles is photoion- 
ized by the obscured active nucleus. We suggest that this type of structure might be 
related to the spatially extended HI absorbers that are detected in front of many high- 
z radio galaxies. We also suggest that we may be witnessing the destruction of the 
extended emission line region. In addition, we investigate the nature of the emission 
line brightness asymmetry in MRC 0406-244: we conclude that this asymmetry is due 
to an asymmetry in the mass of ionized gas, confirming the scenario of McCarthy, van 
Breugel & Kapahi, and leading us to reject that of Copal-Krishna & Wiita. 

Key words: galaxies: active; galaxies: ISM; galaxies: high-redshift; galaxies: jets; 
galaxies: evolution 



1 INTRODUCTION 

Powerful radio galaxies at high redshift (z>2: HzRG here- 
inafter) continue to play an important role in cosmologi- 
cal investigations. They provide a means to identify and 
study the most massive galaxies in the early universe (e.g. 
Rottgering et al. 1994; Seymour et al. 2007), and represent 
an opportunity to examine the possible symbiosis between 
the nuclear activity and the host galaxy (e.g. Nesvadba et 
al. 2006). 

The HzRG MRC 0406-244 (z=2.44: McCarthy et al. 
1991) has been the subject of observations in various wave- 
bands. Radio images obtained at ~5 and ~8 GHz (~16 and 
28 GHz in the rest-frame) show a triple radio source which 
has a relatively small projected diameter (82 kpc), and which 
shows a significant rotation measure (Rush et al. 1997; Car- 
illi et al. 1997). The presence of warm ionized gas, with a spa- 
tial extent of at least 66 kpc, has been revealed through spec- 
troscopy and narrow-band imaging (McCarthy et al. 1991; 
Eales & Rawlings 1993; Rush et al. 1997; Pentericci et al. 
2001; Taniguchi et al. 2001; Iwamuro et al. 2003; Nesvadba et 
al. 2008). The spatially integrated fluxes of the UV-optical 
emission lines are collectively well explained by photoion- 
ization by the active nucleus, with roughly solar metallicity, 



and an internal extinction A„ < 0.5 mag (Humphrey et al. 
2008a). In addition, a spatially extended absorption feature 
has been detected in the two-dimensional spectrum of Lya 
(Pentericci et al. 2001). The two-dimensional spectrum of 
the CIV A1549 line is also consistent with having been ab- 
sorbed (Taniguchi et al. 2001). Using Spitzer photometrey, 
Seymour et al. (2007) have estimated the stellar mass of 
MRC 0406-244 to be io n - 36±013 . 

Hubble Space Telescope (HST hereinafter) images ob- 
tained by Rush et al. (1997) and Pentericci et al. (2001) 
show spatially resolved emission, with a complex figure-of- 
eight morphology aligned with the axis of the radio emis- 
sion. A likely scenario is that this figure-of-eight structure is 
a pair of 'superbubbles' resulting from ambient interstellar 
media being swept up by a superwind (McCarthy et al. 1999; 
Taniguchi et al. 2001) or by the radio source (Nesvadba et 
al. 2008). If indeed this is the case, then MRC 0406-244 rep- 
resents an excellent opportunity to study superwinds and 
feedback related phenomena. 

In this letter, we investigate the nature of the superbub- 
bles, using long-slit spectroscopy of the rest-frame optical 
emission of MRC 0406-244. We also revisit the HST image 
of the rest-frame optical emission of Pentericci et al. (2001), 
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Figure 1. Section of the two dimensional ncar-IR spectrum of 
MRC 0406-244, showing the [OIII] AA4959,5007 doublet. The hor- 
izontal axis gives the spatial offset, in arc seconds, relative to the 
peak of the near-IR continuum emission. Wavelength increases to- 
wards the top of this frame. The four emission line knots studied 
in this letter are marked A, B, C and D. 
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Figure 2. HST NICMOS H-band image of MRC 0406-244, ro- 
tated such that the x-axis and the radio axis (PA=128°) are par- 
allel to each other. As discussed in the text, the spatially extended 
emission is dominated by nebular emission. The emission minima 
in the supcrbubble structures are marked with arrows. The axes 
give the spatial offset, in arc seconds, relative to the peak of the 
H-band emission. 



and examine the nature of the side-to-side line brightness 
asymmetry in the extended emission line region. We assume 
a fiat universe with Ho = 71 km s _1 Mpc~\ Q,a = 0.73 and 
fi m = 0.27, resulting in a spatial scale of 8.2 kpc per arcsec 
at the redshift of MRC 0406-244. 



2 OBSERVATIONS AND DATA ANALYSIS 
2.1 Subaru near-IR spectrum 

The long slit spectrum of MRC 0406-244 was obtained on 
15 September 2000, using the OH Airglow Supressor (OHS 



hereinafter) spectrograph (Iwamuro et al. 2001) at the Sub- 
aru telescope on Hawaii. The total exposure was 6000 s, 
with the J and H bands being observed simultaneously. At 
the redshift of MRC 0406-244, these two bands sample a 
rest-frame spectral range of ~3220-3930 A and ~4290-5240 
A, respectively. The slit was l"wide, and had a position an- 
gle of 132°measured North through East, in order to capture 
the bright emission aligned with the radio source (see e.g. 
Rush et al. 1997; Pentericci et al. 2001). The FWHM of the 
seeing disc was ~0.6". The wavelength scale of the resulting 
long slit spectrum is 8.5 A pixel -1 , and the spatial scale is 
0.111"pixel~ 1 . Iwamuro et al. (2003) give a more in-depth 
discussion of these observations and their reduction. The 
two-dimensional spectrum is shown in Figure 1. 

For the absolute flux calibration, we scaled the spec- 
trum such that the total H-band flux measured in a 
2.3" x 1.0"extraction aperture matches the flux measured in 
a 2.3" x 1.0"simulated slit on the H-band image of Iwamuro 
et al. (2003). We expect the uncertainty in our absolute flux 
scale to be ~10 per cent. We adopt the position of the peak 
H-band continuum flux as the spatial zero, and we assume 
that the active nucleus is at this spatial position. 

Four distinct knots of line emission are seen present in 
the two-dimensional spectrum of MRC 0406-244. For each of 
these knots we have extracted a one-dimensional spectrum. 
Aperture A (2.16"NE - 1.05"NE) captures the knot of emis- 
sion to the North East of the nucleus; aperture B (0.72"NE 
- 0.06"NE) corresponds to the line emission near the cen- 
troid of the continuum emission; aperture C (0.17"SW - 
0.94"SW) captures the brighter of the two emission line 
knots on the South West side of the nucleus; aperture D 
(1.28"SW - 2.28"SW) contains the fainter of the two knots 
to the South West. We measured fluxes for the emission lines 
using the IRAF 0NEDSPEC package. The la uncertainties are 
dominated by the uncertainty in the level of the continuum. 
Our measurements are given in Table 1. Note that the low 
spectral resolution of this spectrum (~1400 km s~ ) pre- 
cludes a detailed analysis of the kinematic properties of the 
emission lines. 



2.2 HST near-IR image 

In addition, we make use of a near-IR Hubble Space Tele- 
scope image of MRC 0406-244 (Proposal ID 7498; PI P. Mc- 
Carthy), which was taken using Camera 2 of the Near In- 
frared Camera and Multi-Object Spectrometer (NICMOS 
hereinafter). The F160W filter was used, giving a rest- 
frame wavelength coverage of 4070-5230A. The plate scale 
is 0.076"per pixel. The image is shown in Fig. 2. For further 
details, see Pentericci et al. (2001), where this image was 
first published. 



3 RESULTS AND DISCUSSION 
3.1 Emission line ratios 

Ratios between emission lines are often used to diagnose 
the physical conditions and the ionization mechanism for 
the ionized gas associated with active galaxies (e.g. Bald- 
win, Phillips & Terlevich 1981). In all four knots, we find 
that the optical line ratios are most naturally explained 
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Table 1. Emission line fluxes and ratios in the knots along the radio axis of MRC 0406-244. Columns are as follows: (1) Aperture 
designation; (2) position in arc seconds of the aperture, relative to the peak of the NIR continuum emission, along the slit; (3)-(7) line 
fluxes in units of 10~ 17 erg s _1 cm~ 2 ; (8)-(10) emission line ratios referred to in the text. 



Ap. 


Pos. 


H/3 


[NcV] 


[Oil] 


[Nelll] 


[OIII] 


[NeIII]/[OIII] 


[OIII]/H/3 


[Oil] /[OIII] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


A 


2.16NE-1.05NE 


^3.8 


^4.6 


<12 


<2.3 


24.6±1.5 


^0.1 


^6.1 


^0.5 


B 


0.72NE-0.06NE 


5.3±0.7 


4.5±1.5 


19.0±2.0 


6.6±0.8 


87.3±4.4 


0.076±0.004 


17±2 


0.22±0.02 


C 


0.17SW-0.94SW 


10.9±1.1 


9.7±1.7 


27.8±1.7 


12.4±1.1 


155.2±7.8 


0.080±0.008 


14±2 


0.18±0.01 


D 


1.28SW-2.28SW 


8.1±1.3 


^3.5 


29.5±2.1 


4.0±1.9 


82.2±4.2 


0.05±0.02 


10±2 


0.36±0.03 



with photoionization by a hard continuum. In particular, the 
[NeIII]/[OIII] ratio (Table 1) is consistent with photoion- 
ization and A„=0 mag (also Humphrey et al. 2008a), but 
is substantially lower than predicted by fast shock models 
(i.e., 0.2-1: Dopita & Sutherland 1996). The [OIII] /H/3 ratio 
is also consistent with photoionization by a hard continuum 
(e.g. Robinson et al. 1987). The detection of [NeV] in aper- 
tures B and C also supports this idea. This conclusion is in 
agreement with that of Taniguchi et al. (2001). It is, how- 
ever, important to recognise that the line ratios themselves 
do not allow us to determine the source of the hard ioniz- 
ing continuum (i.e., the AGN, the metagalactic background 
radiation or radiative shocks). 



3.2 HST NICMOS image 

Extensive imaging in a variety of wavebands has revealed 
spectacular filamentary or bubble-like features roughly along 
the radio axis of MRC 0406-244 (Rush et al. 1997; Penter- 
icci et al. 2001; see Figure 2). While imaging alone provides 
insufficient information to elucidate the nature of these mor- 
phological features, spectroscopic data can be used to assess 
what fraction of this spatially extended emission is emitted 
by warm ionized gas. For each of our spectroscopic aper- 
tures, we first calculated the flux expected due to nebular 
continuum within the spectral range of the F160W filter. 
This we have estimated using the H/3 flux and the nebu- 
lar continuum emission coefficients in Aller (1984), assum- 
ing that the gas has a temperature of 15,000 K. We then 
summed the expected nebular continuum flux with the mea- 
sured flux of line emission, and compared this against the 
total flux measured within the spectral range corresponding 
to the wavelength range of the relevant broad-band filter. In 
the spatially extended figure-of-eight structure (apertures A, 
C and D) , we find that the overall H-band flux is dominated 
by nebular emission (i.e. ^90 per cent). 

While Rush et al. (1997) suggested that the figure-of- 
eight structure has a tidal origin, we prefer the interpretation 
of McCarthy et al. (1999) and Taniguchi et al. (2001), i.e., 
this structure is ambient ISM that has been swept up into a 
bubble by a galactic superwind. Our reasoning for this is as 
follows. In the H-band image, the bright emission knots are 
arranged into two ellipses, one on each side of the brightest 
emission peak, which we assume marks the position of the 
galactic nucleus. Such morphological features can be natu- 
rally explained as gas that has been swept up from the ISM 
of MRC 0406-244 by a galactic wind, and which has since 
cooled at the interface between the wind and the external 
medium. Indeed, as remarked by McCarthy et al. (1999), 



3.1 




Figure 3. Diagram to illustrate our simple scenario, in which 
a uniform spherical shell of gas is irradiated from outside by the 
active nucleus. For several points in the shell, we show the relative 
gas column density measured from that point to the active nucleus 
(see text). 



these features are reminiscent of the wind-blown bubbles 
associated with some ultra-luminous infrared galaxies. In 
addition, at the centre of each of the emission ellipses there 
is a region of essentially zero flux. This can be naturally ex- 
plained if the wind has been relatively efficient at sweeping 
up, or at shock-heating, the ambient ISM, leaving little or no 
cool gas in its wake. Outside of the bubbles, a significant flux 
of emission is detected, which can be explained as emission 
from ISM across which the bubble has not yet expanded. 

In §3.1, we analysed ratios between emission lines, con- 
cluding that the gas in the super bubbles is ionized primarily 
by a hard continuum. As mentioned in §3.1, the line ra- 
tios themselves have not allowed us to discriminate between 
three potential sources of hard ionizing continuum emission. 
In this context, the detailed morphology of the super bubble 
can provide a useful test to discriminate between photoion- 
ization by the anisotrpic radiation field of the AGN, and 
ionization by shocks or by the isotropic meta-galactic back- 
ground radiation (MBR). 

Let us suppose that each superbubble in MRC0406-244 
is a spherical bubble of cool gas, which is optically thick to 
ionizing radiation, which has a uniform distribution of ma- 
terial, and which is irradiated by an AGN located at some 
point just outside of the bubble. Some points on the bub- 
ble will be directly illuminated by the AGN; at other points 
on the bubble, the incident radiation will have been filtered 
or absorbed, because the line of sight to the AGN passes 
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through another part of the bubble. In this simple scheme, 
the region of the bubble that views the AGN through the 
highest column density of gas will be at a position angle of 
±45-90° , measured around the circumference of the bubble 
from the point closest to the AGN (see Fig. 3). The nebu- 
lar emission from this shielded region should be fainter than 
from other regions of the bubble. Hence, in the super bub- 
ble structures of MRC 0406-244, we would expect to see flux 
minima in both bubbles at position angles of ±45-90° , mea- 
sured from the point closest to the AGN. Remarkably, such 
minima are indeed present in the bubble structures MRC 
0406-244 (see Fig 2), and this leads us to suggest that pho- 
toionization by the active nucleus is the dominant source of 
ionization for the EELR of this object. These minima can 
also be seen in the HST F675W and F439W optical images 
of MRC 0406-244 (Rush et al. 1997). 

The [OI] A6300 emission line, which falls into the ri- 
band in the case of MRC 0406-244, would provide a useful 
test for our scenario. In the context of emission line regions 
photoionied by an AGN, the [OI] line is formed in zones that 
are shielded from ionizing UV photons (these photons hav- 
ing been quenched), but where there remains a strong flux 
of soft X-ray photons. In our scenario, these conditions per- 
tain to gas associated with the H-band flux minima in the 
superbubbles. Therein, we would expect to measure unusu- 
ally high [OI]/Ha ratios, if indeed our scenario is an accu- 
rate representation of 'reality'. Thus, K-band spectroscopy 
would be useful to test our scenario. In addition, it would 
be interesting to image MRC0406-244 in sub-mm continuum 
and CO line emission, using the next generation of mm- wave 
interferometers such as ALMA: studying the spatial distri- 
bution of the cold ISM should improve our understanding 
of this galaxy. 



3.3 Physical properties 

Having concluded that the bubble structures are photoion- 
ized by the active nucleus, it is now possible to obtain an 
estimate for their gas density n using 

- Q m 
71 ~ 47rr 2 cU 1 ' 

where Q is the ionizing photon luminosity of the active nu- 
cleus, r is the distance of the gas from the source, and U is 
the ionization parameter. Assuming that the ionizing beam 
of the active nucleus has an opening angle of ~90° (Barthel 
ct al. 1989), the global H/3 luminosity of MRC 0406-244 of 
L/j ( 3=1.6xl0 43 erg s" 1 (Iwamuro et al. 2003; Humphrey et 
al. 2008a) then implies a minimum Q of 1.3xl0 55 s _1 . For 
knot D, at a projected distance of r=12.8 kpc from the AGN, 
and which has an ionization parameter U=0.004 (§3.1), we 
thus obtain a gas density of n^6 cm -3 . For knot C, we 
obtain n^27 cm -3 . These values are consistent with those 
measured by Nesvadba et al. (2008). If we assume that the 
clouds in the bubbles are ionization bounded, then the col- 
umn density of ionized Hydrogen is 

N H+ = Uca recB ^ = 6 x 10 20 cm" 2 (2) 

where a rec _g is the total Hydrogen recombination coeffi- 
cient for case B. In standard AGN photoionization models 
Nhi/Nh > 10" 6 , and hence N H i > 6 x 10 14 cm" 2 . The 



total mass of ionized Hydrogen in the superbubbles can be 
estimated using 

M H+ = u Lh ^ (3) 

from which we estimate Mh+ ~ 6 x 1O 8 M0, assuming that 
emission from the spatially extended superbubble dominates 
the total flux of H/3 (see Table 1 col. 3) . It is important to re- 
alise that in this calculation we have not considered neutral 
or molecular gases, and thus the true mass of the superbub- 
bles could be significantly higher. (See Nesvadba et al. 2008 
for a discussion of the energetics of the superbubbles.) 

3.4 Line brightness and radio lobe length 
asymmetries 

There is a trend amongst powerful double radio galaxies in 
which the spatially extended ionized gas is usually brighter 
on the side of the nucleus which has the shorter of the two 
radio lobes (McCarthy, van Breugel & Kapahi 1991). The 
widely accepted explanation for this trend is that the inter- 
stellar medium is asymmetrically distributed, resulting in 
slower expansion of the radio source and brighter line emis- 
sion on the side of the nucleus where there is more ISM 
(McCarthy, van Breugel & Kapahi 1991). A possible alter- 
native scenario uses motion of the radio galaxy through the 
inter-galactic medium to produce a lobe-length asymmetry 
and an asymmetrical extinction of the optical line emission, 
the latter resulting in a line brightness asymmetry (Gopal- 
Krishna & Wiita 1996). To date, these two scenarios have 
been subjected to few observational constraints, and very 
little new data has been applied to understanding the na- 
ture of these correlated asymmetries (e.g. Ishwara-Chandra 
et al. 1998; Humphrey et al 2007). In this section, we use 
results from our Subaru long-slit spectroscopy to obtain in- 
formation about the nature of the line brightness asymmetry 
in MRC 0406-244. 

As expected, MRC 0406-244 follows the trend, i.e., the 
optical line emission is brighter on the SE side of the nu- 
cleus, which has the shortest radio lobe (e.g. Figs. 1 and 2; 
Rush et al. 1997). We investigate the nature of the bright- 
ness asymmetry in MRC 0406-244, using two apertures in 
which the signal to noise of the emission lines is high, which 
encompass the two brightest emission line knots within the 
slit, and which are approximately equidistant (in projection) 
from the assumed position of the active nucleus: apertures 
B and C (0.72"NE - 0.06"NE and 0.17"SW - 0.94"SW, re- 
spectively; see Table 1). The main results are as follows. 

1. We do not detect any significant difference in the reden- 
ning sensitive line ratio [NeIII]/[OIII] between B and C, sug- 
gesting that there is no significant difference in extinction 
or reddening by dust between knots B and C. Therefore, we 
reject the scenario of Gopal-Krishna & Wiita (1996), which 
attempts to explain line brightness asymmetries as an effect 
of differential extinction by dust. 

2. Moreover, several line ratios that are highly sensitive to 
U, such as [OII]/[OIII] and [OIII]/H/3, show no significant 
difference between B and C. This suggests that there is no 
significant difference between B and C in terms of the local 
gas density iih, or in terms of the local ionizing luminosity 
Q/47rr 2 . (Note: since [OII]/[OIII] is also sensitive to reden- 
ning, this result reinforces result 1.) 
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3. Knot C is a factor of ~2 brighter than knot B: taking 
into account results 1 and 2, the most natural explanation 
is that there is a higher mass of warm ionized gas on the 
SE side of the nucleus than on the NW, in agreement with 
the original scenario suggested by McCarthy, Van Breugel 
& Kapahi (1991), and also in agreement with conclusions 
reached by Humphrey et al. (2007) in the case of the z=2.57 
radio galaxy TXS 0828+193. We note that the radio proper- 
ties of MRC 0406-244 are consistent with this explanation: 
compared to the NW radio hotspot, the SE hotspot has a 
lower radio polarization and higher rotation measure (Carilli 
et al. 1997), consistent with being viewed through a compar- 
atively larger column of magneto-ionized gas. 

Finally, we suggest that high spatial resolution images 
of emission from molecular gas, neutral Hydrogen and dust 
would be useful to explore further the nature of the side-to- 
side asymmetries in this and other powerful radio galaxies. 
For instance, does the cold phase of the interstellar medium 
show a similar brightness asymmetry, and is this correlated 
with the radio lobe-length asymmetry? 

3.5 Implications and speculations 

In MRC 0406-244, the alignment between the superbubbles 
and the biconical ionizing radiation field of the AGN has 
allowed the superbubbles to be seen in emission. It is un- 
clear whether this alignment is by chance or whether the 
radio/optical axis defines the preferred axis of this kind of 
outflow. In this section we ask, what if the superbubbles 
were not within the ionization cones? I.e., what if the bub- 
bles expand beyond, or are misaligned with, the cones, and 
what would be their observational characteristics? 

Being outside of the ionizing radiation field of the active 
nucleus, the superbubbles would likely have a low ionization 
fraction, and thus might not be detectable in emission. How- 
ever, if the covering factor of the superbubble structure is 
sufficient (i.e., >0.1), then the high expected column density 
of the bubbles, estimated to be > 6 x 10 14 cm~ 2 in the case 
of MRC 0406-244 (§3.3), would render them detectable via 
absorption lines such as HI Lya. We would also expect the 
absorption lines to be blushifted relative to the line emission 
from the EELR, due to expansion of the superbubbles, as 
well as spatially extended. 

Such observation properties are similar to those shown 
by many of the HI absorbers that are often detected in front 
of HzRG (e.g. Rottgering et al. 1995; van Ojik et al. 1997; 
Pentericci et al. 2001; Jarvis et al. 2003; Wilman et al. 2004; 
Villar-Martm et al. 2007; Humphrey et al. 2008b), which 
leads us to suggest that the superbubbles in MRC 0406-244 
and the HI absorbers in front of some other HzRG might be 
similar phenomena. 

If the bubbles have swept up most of the cool/warm in- 
terstellar gas, and if the expansion of the bubbles proceeds at 
the velocities implied by the emission line kinematics (sev- 
eral hundred km s~ : e.g., Taniguchi et al. 2001), then the 
majority of the extended gas could leave the ionization cones 
(and host galaxy) on a timescale of ~ 10 8 yr. We speculate 
that we may be witnessing the destruction of the extended 
emission line region of MRC 0406-244, and the galaxy's tran- 
sition to join the 1/3 of luminous radio galaxies that emit 
radio and optical continua, but no detectable emission lines 
(Miley & De Breuck 2008). 
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